Mon. Not. R. Astron. Soc. 000, 1-12 (0000) Printed September 23, 201 1 (MN KTeX style file v2.2) 



Shear-flexion cross-talk in weak-lensing measurements 

M. Viola^'^, P. Melchior^'^, M. Bartelmann^ 

^ Zentrum fur Astronomie, ITA, Universitdt Heidelberg, Alhert-Ueberle-Str.2, 69120 Heidelberg, Germany 

^ Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh, EH93HJ, U.K. 

^ Center for Cosmology and Astro-Particle Physics, The Ohio State University, 191 W. Woodruff Ave., Columbus, Ohio 43210, USA 

"^Department of Physics, Tire Ohio State University, 191 W. Woodruff Ave., Columbus, Ohio 43210, USA 

E-mail: mv@roe.ac.uk 



o 

Oh 

(N 



6 

c3 



> 

O 
(N 

0\ 

cn 

o 



Accepted 2011 September 2 1 . Received 20 1 1 September 1 2 ; in original form 20 1 1 July 20 



ABSTRACT 

Gravitational flexion, caused by derivatives of the gravitational tidal field, is potentially im- 
portant for the analysis of the dark-matter distribution in gravitational lenses, such as galaxy 
clusters or the dark-matter haloes of galaxies. Flexion estimates rely on measurements of 
galaxy-shape distortions with spin-1 and spin-3 symmetry. We show in this paper that and 
how such distortions are generally caused not only by the flexion itself, but also by coupling 
terms of the form (shear x flexion), which have hitherto been neglected. Similar coupling 
terms occur between intrinsic galaxy ellipticities and the flexion. We show, by means of nu- 
merical tests, that neglecting these terms can introduce biases of up to 85% on the F flexion 
and 150% on the G flexion for galaxies with an intrinsic ellipticity dispersion of a^ ~ 0.3. 
In general, this bias depends on the strength of the lensing fields, the ellipticity dispersion, 
and the concentration of the lensed galaxies. We derive a new set of equations relating the 
measured spin-1 and spin-3 distortions to the lensing fields up to first order in the shear, the 
flexion, the product of shear and flexion, and the morphological properties of the galaxy sam- 
ple. We show that this new description is accurate with a bias ^ 7% (spin-1 distortion) and 
^ 3% (spin-3 distortion) even close to points where the flexion approach breaks down due 
to merging of multiple images. We propose an explanation why a spin-3 signal could not be 
measured yet and comment on the potential difficulties in using a model-fitting approach to 
measure the flexion signal. 
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1 INTRODUCTION 

In the last years weak gravitational lensing became a standard tech- 
nique to infer the mass distribution of galaxy clusters. In contrast 
to X-ray or galaxy dynamics, which typically assume hydrostatic 
or virial equilibrium and spherical mass distributions, gravitational 
lensing allows to directly probe the projected gravitational poten- 
tial: Light from distant galaxies gets deflected by the cluster po- 
tential and the galaxy ellipticities are modified by the gravitational 
tidal field of the intervening matter. Hence measurements of the 
galaxy ellipticities allow to directly trace the curvature of the grav- 
itational potential. 

Beyond that, higher-order distortions of background galaxies, 
like curvature or distortions with three-fold rotational symmetry, 
can be related to higher-order derivatives of the gravitational po- 
tential, thereby providing information on the variation of the po- 
tential at smaller scales. By expanding the lens equation to sec- 
ond order, one can identify two combinations of third-order deriva- 
tives of the potential, called F and G flexions. These two fields 
can be related to a skewness in the light distribution and distortions 



with three-fold rotational symmetry in the shape of galaxies, re- 
spectively (Goldberg & Bacon 2005; Bacon et al. 2006; Irwin et al. 
2007). Reliable measurements of these higher-order distortions can 
increase the accuracy with which the inner profile of galaxy clus- 
ters can be constrained or can support the detection of substructures 
(Bacon et al. 2006), whose presence in galaxy clusters is a firm pre- 
diction of the Cold-Dark Matter cosmogony (White & Rees 1978). 
Moreover, at galactic scales the incorporation of flexion constraints 
significantly improves the recovery of halo ellipticity in compari- 
son to shear-only measurements (Hawken & Bridle 2009). 

Two main approaches to measure these higher-order distor- 
tions in the galaxy shape have been proposed so far; shapelets 
(Goldberg & Bacon 2005; Velander et al. 2011) and HOLICs 
(Okura et al. 2007, 2008). Both assume that the measured distor- 
tions can simply be related to the F- and G-flexions, respectively. 
This is true only if the shear can be considered negligible and galax- 
ies are intrinsically circular. We will show that a considerable error 
on the inferred lens power is introduced if any of these assumptions 
is violated. Then, coupling terms of the form (shear x flexion) or 
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(ellipticity x flexion) create additional distortions, which are in- 
distinguishable from the pure flexion-induced distortions and thus 
interfere with the interpretation of the lensing signal. 

In section 3 we derive a new set of equations, relating the 
measured distortion to the lensing fields, including these coupling 
terms. In section 4 we discuss the physical meaning of the newly 
identified terms. In section 5 we derive a simplified set of equations, 
applicable to lenses with non-vanishing shear, which also consider 
the responsivity of the higher-order distortions to shear. These new 
equations are more general than the equation originally presented 
by Okura et al. (2008) and constitute a practical, yet accurate ap- 
proximation to the equations presented by Schneider & Er (2008). 

By means of numerical tests, we show in section 6 that our 
equations are sufficient to reliably describe the lensing-induced 
change of galaxy shapes up to the point where the very concept 
of flexion estimates from individual galaxies breaks down due to 
merging of multiple images of the same source (Schneider & Er 
2008). In section 7 we propose an explanation why the G-flexion- 
related spin-3 field has not been measured yet (Goldberg & Bacon 
2005; Leonard et al. 2007; Okura et al. 2008; Leonard et al. 201 1) 
despite the fact that the G-flexion extends to larger cluster-centric 
distances than the F-flexion, and we comment of the conceptual 
difficulties of using a model-fitting approach to measure the flexion 
signal. We conclude in section 8. 



2 FLEXION FORMALISM 

This section summarises the basic weak-lensing concepts used 
throughout this work. For a complete overview we refer to Bartel- 
mann & Schneider (2001). An isolated lens with surface mass den- 
sity S(0) has the lensing potential 

AG DiD. 



*(0) ^ £J ^hB± f d^e'^{e') \n\e-e'\ 



(1) 



where G denotes the gravitational constant, c the speed of light, and 
Di^s,is are the angular-diameter distances between the observer and 
the lens, the observer and the source, and the lens and the source, 
respectively. In the approximation of a single lens plane, ^ light rays 
are deflected by the angle 



a{e) = V*(0) 



(2) 



which relates the angular positions on the source plane /3 and those 
on the apparent image plane G by the lens equation 

/3 = e-cx{e). (3) 

If we assume that the source is much smaller than the scale on 
which properties of the lens change, we can locally expand the lens 
equation: 

1 
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In many weak lensing applications the expansion of the lens equa- 
tion can be safely truncated at first order, such that the image dis- 
tortion is described by the Jacobian matrix 

1 — K - 7i -72 

-72 1 — K + 71 



A: 



(5) 



^ This approximation provides a very accurate description of most lensing 
phenomena, including flexion (Schafer et al. 20 11). 



where k is the convergence and the 7^ are the two components of 
the shear. If, however, the sources are large or the lens properties 
change rapidly on small scales (e.g. towards the core of a galaxy 
cluster), second-order terms become important." 

At this point, it is convenient to introduce a complex notation 
for the relevant lensing quantities. We start defining a complex gra- 
dient operator (Newman & Penrose 1962; Bacon et al. 2006) 



a = — +i— = 9i +id2, 

Wi 6^2 

which in polar coordinates has the form 
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This representation shows clearly that when d is applied to a spin-s 
quantity,'^ it raises its spin by one. Analogously, d* lowers the spin 
by one. Applying this operator one time to the deflection potential 
(spin-0), we generate the deflection angle 



a = a*, 



(8) 



a spin-1 field. Applying it twice, we generate the spin-0 conver- 
gence field 



^aa** = ^(*,ii + *.22) 



or the spin-2 shear field 
1 



1- 
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2 2' 

while applying it three times generates the so-called T-flexion 
(spin-1) 

j^= ^aaa** = ^[(*,iii + *,i22) + i(*,ii2 + *,222)] (ii) 

or the Q-flexion (spin-3) 

g = ^aaa* = ^[(*,iii-3*,i22) + j(3*,ii2-*,222)] (12) 

fields. We can now use these new fields to rewrite Equation 4 in 
complex notation: 



/3~(1-A.) 
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where the reduced shear and flexions are defined as 



7 



F = 
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, G = 



(13) 
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1 — K 1 — K 1 — K 

This equation relates source plane and observed positions in terms 
of derivative of the lensing potential with uniquely defined spin 
properties. One can show (Schneider & Er 2008) that the Jacobian 
of this transformation is 



dety4~(l-K)^ 



1 



99 



F* 



F+^\9*G + gF* 



2 

+ o{e^ 



9*F + gG' 



(15) 



^ Eventually, even the second-order terms cannot fully describe all lensing 
effects, for instance the occun'ence of giant arcs. In these cases the full 
Equation 3 has to be studied. 

^ We say that a quantity has spin s if it is invariant under a rotation of the 
Cartesian coordinate frame by a rotation angle <f> = 2-k/s and s £ [1,2,...]. 
Vectors are thus spin-1 quantities. 
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3 MEASURING FLEXION FROM MOMENTS 

Gravitational flexion is related to the third derivative of the lensing 
potential (third-order lensing effect) and it manifests itself in the 
skewed and arclike appearance of lensed galaxies. A possible de- 
scription of these higher-order distortions can be given in terms of 
moments of the galaxy surface brightness, 



^ij...k 



<fei{e)e,ej 



(16) 



where So denotes the flux (0-th moment). For instance, a common 
ellipticity estimator can be formed from second-order moments, 



2i<^ 



+ 62 



(Qll-g22) + 2^Ql 
Qii + Q22 



(17) 



where a and b are the semi-major and semi-minor axes of the el- 
lipse. Consequently, it is possible to form combinations of third- 
order moments such that they exhibit spin-1 or spin-3 properties, 
which we then seek to relate to the flexion fields. In this work, we 
follow the notation introduced by Okura et al. (2008) to describe 
these higher-order distortions: 



5^^J' cFei{e)e' = 



, ^ (Qui + Q122) + i{Qii2 + Q222) 
(Qui — 3Q122) + i(3Qii2 — Q222) 



(18) 



f = 



C - 23C* - 3*<5 - 2F*7? 



[1 - 43^(5*7,)] (1-^)^ 

'l^^ ^Gr,* - ^G'\ + 3F*g + ^Fg*r, + ^{F*g*\) 

+ IcgX* + ^-Fgr,* + j4F*x + \Fg*X + 9G*X + 3F 



+ Gg* - 4F*3 - ^gFx* + \gx* - 2>F^{gx* 



(20a) 



and 

5^ = 



S-SgC- ^F-q- \f*\~ ^G 
a-, 2 ' 4 4 



[l-4K(»*r,)](l-^)^ 

+ &Fg - ]^Fg*\ + 4F*ot + ^-Gg^* - \g* g\ - ^Gg^v 

+ sJIfx+ \Gg*x - \F*gx - 3Fg + ]^Ggx' 

' (20b) 

where r\ and A are the following dimensionless spin-2 and spin-4 
quantities, based on 4th-order moments of the light distribution. 



(21) 



(Qiiii — Q2222) + 2i(Qiii2 + Q12 



These quantities have spin 1 and spin 3, respectively. The normali- 
sation factor has spin and it is defined as a combination of fourth- 
order moments. 



_ (Qiiii — 6Q1122 + < 



+ 4i(Qiii2 - 



(22) 



C= J- I d^ei{6){66*f 



)llll+2Qll22+Q2222 , (19) 



and 



Okura et al. (2008) derived the lensing transformations of ( and <5, 
considering only terms of first order in the flexion fields. Specifi- 
cally, terms of the kind shear x flexion were discarded, which is 
equivalent to assuming a vanishing shear As we shall show in sec- 
tion 6 using the example of a realistic NFW halo, flexion will nec- 
essarily be measured in regions of the lens where the shear is not 
negligible (cf. Figure 2). Schneider &. Er (2008) derived a complete 
set of transformations of 3rd- and 4th-order moments under lensing 
in a somewhat different notation, but noticed biases of their flexion 
estimates that remained unexplained. 

We therefore re-derive the transformations of C^ and 5 under 
lensing, keeping all terms of first order in the shear, the flexion, and 
the product of shear and flexion. The calculation proceeds along 
the following scheme (and can be reviewed by the interested and 
courageous reader in Appendix A): 

(i) Compute the transformation of 3rd- and 4th-order moments 
under lensing according to Equation 13. This computation is done 
neglecting the centroid shift induced by lensing. 

(ii) Correct for the lensing induced centroid shift. 

The results of these calculations show that at first order in the flex- 
ion and shear fields the transformations of 1^ and 5 under lensing 



/i; 



Tr{Qf _ (Qii + ( 



e 
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(23) 



is a dimensionless spin-0 quantity that describes the radial profile 
of the source. The quantities on the left-hand side, ^'^ and 5^, de- 
note the intrinsic, pre-lensing spin-1 and spin-3 distortions of the 
source. The quantities on the right-hand side, C,, S, rj. A, x, denote 
the observable distortions. Hence, Equation 20 relates the intrinsic 
distortions with apparent, lensing-induced ones {(, 5) via differ- 
ent combinations of the lensing fields, involving moments up to 
order four of the light distribution. These equations are conceptu- 
ally equivalent to the more familiar one for spin-2 distortion re- 
lating the intrinsic galaxy ellipticity to the shear and the observed 
ellipticity (e.g. Seitz & Schneider 1997). Note that we dropped all 
moments higher than fourth order in the derivation of Equation 20 
as we regard them to be practically immeasurable. Even then, the 
occurrence of terms like 77 and A is worrying given the limited sig- 
nificance of typical background galaxies. For this reason, Leonard 
et al. (2007); Okura et al. (2008); Leonard et al. (2011) have pro- 
posed simplified variants of these equations, exploiting the follow- 
ing additional approximations: 

• ri and A are small and thus negligible; 

• The shear and the convergence are small. 

Incorporating these simplifications, a linear relation between ^ and 
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5 and the flexion fields can be found, which allows us to directly 
solve for 






(24a) 



(24b) 



We shall investigate the validity of this simplified formula in sec- 
tion 6. 



4 MOIfPHOLOGICAL DEPENDENCE OF FLEXION 
MEASUREMENTS 

In this section we seek insight into the physical meaning and im- 
portance of the terms appearing in Equation 20. This will guide 
our approach to simplify these equations such that we can come up 
with practically useful, yet accurate flexion estimators. 



4.1 Galaxy ellipticities x and ?? 

In Equation 20, the ellipticity of a galaxy is described by two quan- 
tites, r] and x- Seitz & Schneider (1995) showed that x transforms 
under lensing as 



X 



X~2g + g\* 



(25) 



l + lgP-2K(3X*) 
It is easy to see that at first order in the shear 

(X) ^ 25(1 - (xx)) , (26) 

where the term in brackets is the so-called responsivity, and (xx) 
denotes the variance of the x-distribution. Even though x is not an 
unbiased estimator of the shear ~ hence the need for the first-order 
responsivity correction factor - it is normally used in weak-lensing 
applications since it has simpler noise properties compared to the 
unbiased shear estimator 

a — b 



2i<p 



~ a + b 
which is related to x according to 

2e 



(27) 



(28) 



The other quantity describing the ellipticity of the object is rj. 
It is a spin-2 quantity, defined by a ratio of fourth-order moments of 
the light distribution. Assuming an object with elliptical isophotes, 
we show in Appendix B that rj is related to e by 



V 



36(1 + i 



2icj> 



l + 4eP 



(29) 



We also compute how rj transforms under the lens equation and 
derive its responsivity analogously to the responsivity of x, 



(v) -39\ 1- oivv) 



(30) 



where {rjrj} is the variance of the 77-distribution. Note that Okura 
& Futamase (2009) introduced rj as higher-order shear estimator 
(Equation 12 in their paper), but did not consider its responsivity. 

To quantify the responsivity corrections for x and rj we need 
to compute the dispersions of the respective distributions. This is, 
however, hampered by their non-linear relation to the proper el- 
lipticity e. Even if we follow the usual assumption that both com- 
ponents of e have Gaussian distributions, this property is lost for 




Figure 1. Distributions of Xi (green) and rji (blue) under the assumption 
that Ei (red) has a Gaussian distribution with standard deviation CTe = 0.3. 
The relations between x, rj and e are given in Equation 28 and Equation 29. 



Table 1. Values of the ji parameter, as defined in Equation 23 for three 
different Sersic index and three different ellipticities. ris = 0.5 correspond 
to a Gaussian profile, while rij, = 4.0 corresponds to a de Vaucouleurs 
profile. 



|e| = |e|=0.3 |e| = 0.6 



=0.5 0.5 

=2.0 0.1 

=4.0 0.005 



0.43 
0.08 
0.004 



0.35 
0.06 
0.002 



Xi and rji. In Figure 1 we show what these distributions actually 
look like under the assumption that e^ has a Gaussian distribution 
with (Te- = 0.3. Calculating the analytic form of the X" and re- 
distributions and their variances is thus not trivial, a numerical in- 
tegration on the other hand provides us with sufficiently accurate 
fitting formulae, valid in the range (Ttj G [0.2,0.4]: 



and 



(xx) = -l-87a? + 2.04O-. + 0.02, 



(Tjrj) = -2.04(7, + 1.96a, + 0.13. 



(31) 



(32) 



4.2 Galaxy concentration /i 

The quantity ji defined in Equation 23 is the ratio of two quantities 
with spin-0, Tr{Q)^ and 5, both describing the size of the object, 
one in terms of 2nd-order moments, the other in terms of 4th-order 
moments. This ratio therefore encodes the steepness of the object's 
surface brightness profile. In other words, ji can be seen as a con- 
centration parameter: the smaller ji, the more concentrated the ob- 
ject is. The intrinsic value of /j, depends on the morphology of the 
galaxy which is used to measure the flexion. It is easy to show an- 
alytically that in the case of a top-hat surface brightness profile, ji 
equals | . A numerical integration of more realistic radial profiles 
(e.g. Sersic profiles) reveals that ji typically falls into the range 
[0.001, 0.5] (cf. Table 1). 

In general /i depends not only on the galaxy's Sersic index but 
also on its intrinsic ellipticity. For a fixed Sersic index, ji decreases 
with increasing ellipticity. The /i-terms in Equation 20 arise from 
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the flexion-induced shift of the centroid position (cf. Appendix A). 
More concentrated galaxies are less susceptible to this shift, hence 
the observable flexion effects are stronger for those. Flexion effects 
are therefore most easily observable for circular galaxies, but also 
for such with low ^l. 



4.3 Spin-4 distortion A 

The remaining quantity A has spin 4 and is thus related to a dis- 
tortion of the galaxy with fourfold rotational symmetry. We expect 
the pre-lensing value of A to be small given the fact that square- 
shaped galaxies appear to be less numerous in Nature (Kayser et al. 
1989). Moreover, its lensing induced change can only be related to 
terms like GF, g^ , gGF* or any other product of lensing fields 
corresponding to a spin-4 field. This implies that all terms of Equa- 
tion 20 containing A are of higher order in the lensing fields and are 
therefore neglected in what follows. 



5 FLEXION ESTIMATORS FOR STRONGER LENSES 

As we argued above, the assumptions made at the end of section 3, 
which led to the linear flexion estimators given in Equation 24, 
appear to be too restrictive, especially for cluster-lensing applica- 
tions. In particular, shear and convergence cannot be neglected at 
distances of [0.1 — 1] Mpc//i from the cluster centre, where flexion 
has been measured so far (Leonard et al. 2007; Okura et al. 2008). 
There the higher-order terms we kept in the derivation of Equa- 
tion 20 become important. At larger scales, where the assumptions 
of small shear and convergence hold, the signal is very weak due to 
the rapid decrease of the strengh of flexion as a function of scale, 
and therefore very difficult to measure. So far, no detection of a 
flexion signal has been claimed on scales larger than 1.5 Mpc h^^. 
To first order in the shear, the main corrections to the simple flexion 
estimator given in Equation 24 come from terms of the form (shear 
X flexion). 

While the change of spin-1 or spin-3 distortions under lensing 
given in Equation 20 is valid for any single object, a direct solution 
for the lensing fields is impossible since we do not know the intrin- 
sic values (^° and 5^} However, if these distortions obey isotropy, 
they should average to zero as long as the background source en- 
semble is not affected by lensing. This means that the left-hand 
sides of Equation 20a and Equation 20b should vanish in the mean. 
We proceed to calculate flexion estimators by averaging the right- 
hand sides of Equation 20a and Equation 20b - assuming that the 
lensing fields are constant for all galaxies considered in the average 
- according to the following scheme: 

(i) Taylor-expand the denominator to first order in g; 

(ii) Average all terms in the equation obtained this way; 

(iii) Drop all terms of second order in g; 

(hi) Use Equation 26 and Equation 30 to link (x) and {jf) to g\ 

(v) Drop all the terms containing products between g and (x) or 
{77) appears since they effectively scale as g^; 

(vi) Neglect all A-terms as they are of at least second order in 
the lensing fields; 



(vii) Express {(^) in terms of F using Equation 24a and (5) in 
terms of G using Equation 24b. This approximation is sufficient 
since any further correction involves terms of second order in g\ 

(viii) Neglect terms like {^rfR,{gx*)), or in general terms like 
{abc) , where a, b and c represent some morphological properties of 
the galaxy. These terms appear because of the Taylor expansion of 
the denominators in Equation 20a and Equation 20b. Each of these 
terms can give a correction to the equations above at the per-cent 
level, but practically their measurement would be extremly noisy. 

Following these steps, we can reduce the complicated form of 
Equation 20 to a much more compact one: 



(C) 



\-^^^)-\{^^\x^)\F 



\-\{^^\x?)-{^^)\Gg' 



+ 



15 



3(m1x1 )-2W \F*g 



{5) ^\{l- (mIx1'»G+ f f - |(m1x|') \9F 



(33a) 



(33b) 



For a given set of lensing parameters F, G and g, the terms on 
the right-hand side now only depend on the mean concentration and 
the ellipticity dispersion of the galaxies used to measured flexion. 

Even at linear order in the flexion fields (i.e. if the shear- 
flexion cross terms are neglected), the ellipticity dispersion of the 
galaxies affects the flexion estimators. This can be seen looking at 
the first term of the right-hand side of the equations above. The 
reason is that the lensing-induced centroid shift is proportional not 
only to the lensing fields but also to the ellipticity of the object: 
the higher the ellipticity, the more the centroid is shifted. The cor- 
rection terms we computed, which do not appear in any previous 
work, take this effect into account. For a Sersic-type galaxy with 
Sersic index Us = 1.5, we found numerically (/i|xP) — 0.05. 

It is interesting to note how galaxy morphology affects the 
strength of the measured distortion. In the extreme example of cir- 
cular top-hat galaxies {p — 0.75), {Q becomes almost completely 
insensitive to F-flexion, but remains susceptible to Gg* and F*g. 

Equation 33 does not provide flexion estimators, but describes 
the average measured spin-1 and spin-3 distortions in terms of lens- 
ing fields and galaxy properties. In principle, assuming that the 
shear is known - e.g. from measurements of x - we could invert 
Equation 33 in order to obtain F and G flexion estimators in terms 
of (", S, ^ and x'- This inversion, however, is in general not neces- 
sary since the comparison with a theoretical model can always be 
done at the level of measured distortions. For this reason we shall 
always refer to distortion estimators in the following rather than 
flexion estimators. 



■* For the same reason it is impossible to infer the shear from a measure- 
ment of the ellipticity of a single galaxy since its intrinsic ellipticity is un- 
known. 



"" Near cluster centres, the estimation of g from galaxy ellipiticities is prob- 
lematic in itself because shear measurements are affected by the presence 
of flexion (Schneider & Er 2008). 
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r [Mpc//i] 



Figure 2. Reduced shear g and reduced flexions G and F (top to bottom 
lines) used in the simulation. The cluster follows an analytic NFW-profile, 
with masses of 10^^ A/q (top panel) and 10^* Mq (bottom panel) at red- 
shift z = 0.5. Background galaxies are located at Zg = 2. 



6 TESTS 

In this section, we assess the performance of the new distortion es- 
timators presented in the previous section by means of dedicated 
simulations. The galaxies in our simulations follow an elliptical 
Sersic-type profile (Sersic 1963), 



I{r) = 7o exp 



l/n 



(34) 



where Re denotes the radius containing half of the flux, Wg the 
Sersic index, and b^^ an integration constant that depends on n^ 
(cf. Graham & Driver 2005). This type of profile is identical to a 
Gaussian for Us — 0.5 and is steeper in the centre for ris > 0.5. 

For our test we assumed Us = 1.5, which represents the aver- 
age value for rather bright galaxies in the COSMOS field (Sargent 
et al. 2007) and we chose two different radii. Re — 2 pixel and 
Re = 8 pixel to investigate the size-dependence of flexion estima- 
tors. We fixed the pixel-size to be 0.2 arcsec. We lensed the galaxy 
according to Equation 13, rotating the reference frame such that 
the imaginary component of the shear vanishes. The strength of 
the lensing fields was drawn from analytic NFW profiles (Navarro 
et al. 1997; Bartelmann 1996; Bacon et al. 2006) with mass lO" 
or 10^^ Mq and concentration c — 4.67 and c — 3.66, respec- 
tively (Neto et al. 2007), both at redshift z — 0.5 with sources at 
Zs = 2.0. The reduced shear and flexions for this density profile 
are shown in Figure 2. The crucial point of this setup is that all the 
relevant lensing fields are obtained from a realistic cluster profile, 
in particular none of them vanishes in the investigated region. The 
range of field strengths for the two halos is specified in Table 2. 

Before the moments were measured, we recentered the object 



Table 2. Minimum and maximum strength of the reduced shear and the 
reduced flexion used in the simulation. These values are drawn from the 
analytic NFW profiles of Figure 2 in the range r £ [0.15 — 2]Mpc//i. 



M = lO^Hf, 







M = 10"M, 
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Reduced shear g 
Reduced F [arcsec^ 
Reduced G [arcsec" 



0.086-0.30 
0.00065-0.027 
0.0013-0.029 



0.026-0.11 
0.0001-0.020 
0.0005-0.023 



to account for the lensing-induced centroid shift'' . We measured the 
moments of the object and computed S and (, as well as /i and the 
ellipticity x- In this way we could describe the averaged measured 
distortions ((") and (S) in terms of the lensing field and the morpho- 
logical properties of the galaxy sample. We then computed the ra- 
tio of the actually measured distortions to the estimated distortions 
for the two halos, obtained either from Equation 24 (linear estima- 
tor) or from Equation 33 (new estimator derived in this paper). It is 
worth noting here that our simulated galaxy images are not affected 
by noise or convolution with the Point-spread Function of the tele- 
scope, thus we can directly measure unweighted moments of the 
post-lensing galaxy shape. 

In Figure 3 and Figure 4, we show results for circular and 
realistically elliptical galaxies, respectively. In both cases, zero in- 
trinsic flexion was assumed. Both figures show that the linear dis- 
tortion estimators can capture only part of the measured signal, es- 
pecially at small scales or in general where shear and flexion are 
stronger. In particular, the difference between the linear estimator 
and the measured signal can be up to a factor of 3.7 for (5) or 
2.1 for {() for large value of the lensing fields (i.e. very close to 
the centre of a massive cluster) and intrinsically circular sources. 
If ellipticity dispersion is included in the simulation, the measured 
signal is up to 2.5 (S) or 1.7 (Q times larger than what the linear 
estimator predicts. The reason for this substantial underestimation 
of the measured signal by the linear estimators is that C, and 5 are 
just a measurement of a spin-1 and a spin-3 field, which can be 
generated by whatever combination of lensing fields and galaxy 
morphology measures whose product has either spin-1 or spin-3. 
Hence, Equation 33a and Equation 33b, incorporating the coupling 
between shear/ellipticity and flexion, provide a much better de- 
scription of the measured signal, especially if the shear is not very 
small or the ellipticity dispersion does not vanish. The linear esti- 
mator, which generally predicts lower distortions, performs better 
for intrinsically elliptical objects because lensing-induced distor- 
tions are smaller for elliptical rather than circular objects. 

However, we notice a residual percent-level bias in our new 
estimators, which we can explain with the approximations made 
in deriving Equation 33: Since this bias seems not to depend on 
the strength of the lensing fields, we must conclude that it origi- 
nates from neglecting terms accounting for the morphology of the 
galaxy sample over which the average is taken (cf. item (viii) of our 
derivation scheme). As more such combinations can be found for 
spin-1 than for spin-3, the estimator for (" is more strongly affected 
than the estimator for 5. 

The general conclusion from our tests is that the description of 
spin- 1 and spin-3 distortions in terms of shear, flexion and galaxy 
morphology that we derived in Equation 33a and Equation 33b 



^ This is a crucial step since the lensing-induced centroid shift has been 
derived assuming that the centroid is at the coordinate origin of the lens 
plane. 
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Figure 3. Ratio between the measured spin-3 signal (5) (upper panels) and the spin-1 signal {C,) (bottom panels), the linear estimator (black) described in 
Equation 24 and the new estimator presented in Equation 33a (blue) and Equation 33b (red). The effective radius R^ of the galaxies is set to 2 pixels (solid 
lines) or 8 pixels (dashed lines). All galaxies were assumed to be intrinsically circular The two dotted vertical lines show the scale radius where spatially 
connected multiple images occur (dependent on the galaxy size Re). The left panels refer to the halo of M = 10^* Mq while the right panels to the 
M = 10^^ Mq halo. The value of the reduced shear is indicated by the top horizontal axis. 



is valid throughout a wide range of cluster-centric distances. Us- 
ing them we obtain essentially unbiased estimates of the lensing- 
induced distortions up to the point, where the very concept of flex- 
ion breaks down: when multiple images of an extended source 
merge (Schneider & Er 2008). Such a case is exemplified in Fig- 
ure 5 and the appropriate hmits are indicated by dotted vertical lines 
in Figure 3 and Figure 4. Beyond this limit, any flexion measure- 
ment is meaningless. 

Finally, the bias is also a function of the galaxy size. Since 
flexion is a dimensional quantity with unit length"^, the accuracy 
of the distortion estimates does not only depend on the strength 
of the lensing fields but also on the galaxy size: Larger galaxies 
react more strongly to flexion. Hence, for large objects, the higher- 
order corrections to Equation 33 neglected in this paper become 
important for lower values of the input flexion field, or farther away 
from the cluster centre. Also, merging multiple images then occur 
at larger distances. These effects can be seen in Figure 3, where we 
show results from sources with Re = 2 and Re = 8 pixels. 



7 IMPLICATIONS 

Despite several attempts and opposed to the spin-1 flexion distor- 
tion, the lensing-induced spin-3 field has not been measured so far 
(Goldberg & Bacon 2005; Leonard et al. 2007; Okura et al. 2008; 
Leonard et al. 201 1). From a theoretical point of view, this comes as 
a surprise because the G-flexion field generated by a typical NFW 
halo is stronger than the F-flexion field (cf. Figure 2). However, 
when we compare the strength of the measured {() and (5) signals 



in Figure 6, we see that in the radial range r £ [0.2, 1.0] Mpc/h, 
the estimator ((") is between 20 % and 50 % larger than (S), hence 
opposite to the behavior of the F and G fields. 

Furthermore, assuming that the noise on the measured mo- 
ments is Gaussian, we expect the ratio between signal-to-noise ra- 
tios of (^ and 5 to scale as 

where Si denotes the strength of the measured spin-1 (Q and spin- 
3 (S) field, as shown in Figure 6. The pre-factor originates from 
the definitions of (^ and S in Equation 18. This implies that for a 
given noise level, the S/N for i^ can be up to 4 times the one for 5. 
This simple argument could be one explanation why no spin-3 field 
around otherwise prominent lenses has been detected yet. 

On the other hand, the scale where the lensing fields have a 
significant amplitude does not only depend on properties of the 
lens, but also on the dispersion of the intrinsic distortions ^'^ and 
5^ (as well as e°). Since three-fold rotational symmetry is less com- 
mon in galaxies than skewness, we expect as" < ct^s . Moreover, 
as pointed out by Velander et al. (2011), "light leaking" from bright 
objects can also contribute to the measured spin-1 distortion of the 
object but hardly to the measured spin-3 distortion, possibly caus- 
ing a bias in the conversion between ( and the lensing fields. Thus, 
despite the lower S/N, the lower intrinsic dispersion and the dif- 
ficulty to generate spurious spin-3 distortions could render the S 
signal more accessible and cleaner than (^. 

Another interesting point we want to discuss here concerns 
model-fitting techniques for estimating the strength of the flexions. 
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Figure 4. Equivalent to Figure 3 for galaxies with R^ = 2 pixel. The average is taken over a galaxy population with ellipticity dispersion cr^ = 0.3. The error 
bars show the l-cr error of the mean computed by averaging f and S over 1000 galaxies. 




Figure 5. Break-down of the flexion concept. Left panel: Single image of a 
circular source wifli fie = 2pxlensedby3 = (0.18, 0),F = (0.0008,0), 
G = (0.0011,0). Right panel: Multiple images of the same source lensed 
by g = (0.324, 0), F = (0.0067, 0), G = (0.007, 0). 



The standard approach is to model a circular galaxy with a given ra- 
dial profile, to apply spin-1, spin-2, and spin-3 transformations, and 
to compare this distorted model with the actual image in order to 
determine the amplitude of the transformation (e.g. Velander et al. 
201 1). As we have shown here, a spin-1 and a spin-3 distortion can- 
not be simply related to flexion, thus the interpretation of the ampli- 
tude of the best-fitting transformation in terms of lensing fields is 
not straightforward. While shearing an elliptical galaxy is equiva- 
lent to applying a magnification, a rotation and a pure spin-2 distor- 
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Figure 6. Measured (f ) and (5) signal as a function of cluster-centric dis- 
tance for circular galaxies with Re = 2 pixels and the halo with 10^^ Mq. 
The dotted vertical line represent the radial scale where connected multiple 
images appear. 



tion to a circular galaxy (in the limit of small shear: Kuijken 2006), 
the same is not true for the flexions, even if F, G and g are small. 
The reason for this somewhat counterintuitive behavior stems from 
the dependence of the response to flexion on the galaxy ellipticity 
(cf. Equation 33a and Equation 33b). In other words, shearing and 
flexing an elliptical galaxy is in general not the same as applying 
magnification, rotation, spin-1, spin-2 and spin-3 distortions to a 
circular object, just because the amount spin-2 distortion differs in 
the two cases. 

In Figure 7, we compare the strength of the average spin-1 
and spin-3 distortions of intrinsically circular or elliptical galax- 
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Figure 7. Ratio between the average spin-1 and spin-3 distortions for in- 
trinsically circular and elliptical Sersic-type sources. For the circular galax- 
ies, the spin-2 distortions was adjusted such as to generate the entire post- 
lensing ellipticity. Plotted are average and l-cr errors of 1000 galaxies with 
ellipticities sampled from a Gaussian distribution with (Tj = 0.3. 



ies. For the former, the effective shear was chosen such as to fully 
account for the apparent post-lensing elhpticity (thus combining 
intrinsic and actual lensing-induced ellipticity) of the galaxies used 
in the latter case. At large radii, employing a circular source model 
leads to a 40% overestimation of the true spin-1 and spin-3 distor- 
tions. In contrast, at small radii, the actual signal would rather be 
underestimated. The complicated impact of the assumption of cir- 
cular source models on the flexion-distortion estimates would thus 
require the inclusion of the intrinsic source ellipticity into the set 
of model/lensing parameters. Besides increasing the numerical de- 
mands of doing so, the strong degeneracy between intrinsic and 
apparent post-lensing ellipticities, however, impedes such an ap- 
proach without additional assumptions on the strength of the shear. 



8 CONCLUSIONS 

In this work we show that measurements of spin-1 (Q and spin-3 
distortions (5) of galaxy shapes in terms of surface-brightness mo- 
ments cannot be related in a straightforward way to the F- and 
G flexions caused by gravitational lensing, in particular not for 
strongly elliptical galaxies or near the centers of galaxy clusters. 
In the latter case, the reduced shear is not small, thus terms cou- 
pling shear and flexion produce non-negligible spin-1 and spin-3 
distortions, which contribute to the distortions induced by the flex- 
ion alone. A similar effect applies in cases where the observed el- 
lipticity is not small. These coupling terms imply that in general the 
measured signal will typically be larger than the naively expected 
flexion signal. 

Furthermore, lensing causes a centroid shift that is approxi- 
mately linear in the flexion fields and it depends on the ellipticity 
and the concentration of the source galaxy. Hence, when flexion- 
related distortions measured from individual objects are averaged 
over a source galaxy population with given distributions of elliptic- 
ity and concentration, the measured signal will depend not only on 
the strength of the lensing fields but also on average galaxy proper- 
ties. 

Our new derivation of the lensing-induced changes of the spin- 
1 and spin-3 distortion estimators accounts for all of the effects 
mentioned above. This work extends that of Okura et al. (2008) 



- yielding a linear relation between flexion and spin- 1 and spin-3 
distortions (cf. Equation 24) - in two ways: 

• We include shear-flexion coupling terms in the transformation 
of the flexion estimators (^ and 5; 

• We compute how the average of 5 and ( over a population of 
galaxies is affected by the galaxy ellipticity dispersion. 

We provide two more detailed equations relating (i^) and (5) 
with the reduced shear and the reduced flexion. These equations 
take the contribution to the measured signal, originating from con- 
centration and ellipticity dispersions of the galaxy population, ex- 
plicitly into account. Treating the dependence on galaxy morphol- 
ogy is important, since it impacts on the distortion estimates even 
if shear-flexion cross-terms can be neglected. 

By means of several numerical tests, we show that simple lin- 
ear estimators can lead to biases of up to 85% and 150% for (^ and 
5, respectively. This bias depends in general on the strength of the 
lensing field, on the intrinsic ellipticity dispersion and on the con- 
centration of the galaxies. Moreover, we showed that the descrip- 
tion of the spin-1 and spin-3 distortions that we derived in this pa- 
per is valid (with a bias ^ 7% (spin-1 distortion) and ^ 3% (spin-3 
distortion)) even close to the point where the flexion concept breaks 
down due to the merging of multiple images. 

We finally find that the measurement of the spin- 1 distortion 
is always stronger (at scales where flexion is usually measured) 
than the spin-3 distortion. This, combined with the different noise 
properties of the ( and 5 estimators, could be an explanation of why 
a spin-3 field has not been detected so far. 
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APPENDIX A: TRANSFORMATION OF C AND S 



Using the definition of (^ from Equation 18 and the second-order lens equation (Equation 13), we can relate the tensing fields with the spin-1 
distortion on the source plane 



(f9I{e) ( e'-^e" - g*9^ + ^F*g*9'^ - ^G*6l* - 2F*e'^9* + ^Fg*9^9* - 2g9{9*f - 7^^(6161*)^ + 3F*g{99*f + 1:Fg9{9*f 



\9G{9*f 



(Al) 



In the same way we compute how 5 transforms 



5" ^ 



e 



d^9I{e) det ^ ( 61 - g9* - \f*9'^ - ^FeS* - \G{9*f 



^ A <f9I{9) i 9^ - Sg^^e* - ^F9^9* - jf'*6''' - ^G{99*f + 6Fg{99*f - ^^3*6''' + 4F 



e 



^Gff*e^e* I , 



+ lGg{9*f9 (A2) 



where det yl is defined in Equation 15, and 

(," = / d^9I{e) detA[9- g9* - jF*9'^ - ^^FOe* - ^G(6i*)^ 



- /e - If^ f - ^F*e*e - ^G*(e)- 



(1 - kY / d^6i/(6') 



)^ - 2/61^61* 



(A3) 



{l-Kfai-m{gV*)) 



In this derivation, we neglected all moments of higher than fourth order since they are practically immesurable. The combination of these 
two equations leads to Equation 20a and Equation 20b without all the terms containg /i which come from the inclusion of the centroid shift, 
which is calculated in the following. 

Centroid shift 



Normally moments of the light distribution are computed with respect to the centre of light. This means mathematically: 

e= ^ [ d^99I{e) = . (A4) 

JO J 

However, the centroid is not mapped to the origin of the source plane as we assumed in Equation Al and Equation A2. We derive here the 

mapping between the centroid in the source plane and in the lens plane up to first order in tensing fields: 



(1-^)^ 



/.= -yd^/3/3r(/3).A_Z_/,^,,(,),,_^ 



^F*9^ ^ \F9e* - \G{9*f ](i-9(f*^ 9* F + gG* 



9*1 F + 



g*G + gF 



2 (Tr(0))(l - ^) I - ^F*x - \Fg*X - \gG*x ~\f - ^Gg* + \f* g + gFx* - \gx 



(A5) 



Using the form of the centroid shift calculated above we can compute: 



c 



^ey 



and 



^oe 



5* 



d^/3(/?-/3)^(r-/?*)r(/3) 



/d2/?(/3-^)^r(/3) 



5*0^ 



/' 



^7 



/• 



/ d2/3/»(/3)((/3 - ^)(/3* - fi*)y 5(1 - A.)«(l - 4K(ot*) 



d^l3r{l3)(3^l3* - P* / d^pr {13)13^ -213 / d'' pF {13)1313 



d^pF{p)p^ ~?,p / d^pr{p)p'^ 



(A6) 



(A7) 



where we made use of the fact that to the order we are interested in, ^° is unaffected by the centroid shift. The first integral in both expressions 
can be replaced with the expressions computed in Equation Al and Equation A2, while the second and the third integral, including the 
corrections due to the centroid shift, can be replaced by: 



d^pr{P)p^ ~ (1 - ^f I d^9m{9^ - 2g99*) ~ (1 - ^)^{TrQ){x - 2g) 



(A8) 
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<fpr{P)PP* ~ (1 - ^)\TrQ){l - 2K(gx*) 



since all terms proportional to flexion or g^ will be dropped anyhow due to the multiplication with /3. Then we have: 



c = c 



and 



Tr{Q) 
C(1-4K(577*))(1-A.)2 



Tr{Q) 



(X - 2ff)r + 2(1 - 2^{gx*))P 



■i(xP) - &9P 



C(1-4K(ot*))(1-ac)2 
which lead to Equation 20a and Equation 20b, where the upper index "t" has been dropped. 



(A9) 



(AlO) 



(All) 



APPENDIX B: HIGHER-ORDER SOURCE ELLIPTICITY 

We derive here the relation between e ellipticity and t] ellipticity for a source with elliptical isophotes. Without loss of generality we assume 
the source profile to be described by an elliptical gaussian with ellipticity e. Furthermore we select the reference frame such that £2 = 0. We 
can then derive this general relation for the moments of order i + j: 



(i+j) 



dxdyx'y'eyip\ ^l{l ^ efxA e^pi --{1 + efy'^ - ^('t')r('t') 



(l-e)'+i(l + eV+i 



(Bl) 



Using the expression for rj in terms of 4th order moments given in Equation 21, we find 

3e(l + 6^) 
'^ 1 + 4e2 + £4 

Moreover, we are interested in computing the responsivity of rj. For doing so, it suffices to compute the transformation of 77 under lensing 
up to first order in the shear. After some algebra it is possible to show that the following relation holds: 



(B2) 



V 



1-4K(ot*) 



and therefore 



iv) -3g\l- -{rirj) 



(B3) 



(B4) 
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